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Abstract 



q 

' We consider a state estimation problem where observations are made by multi- 

' pie sensors. These observations are communicated over a lossy wireless network to 

a central base station that computes estimates via a Kalman filter. The goal is to 
determine the optimal location of the base station under a certain class of packet 
loss probability models. It is shown in the two sensor case that the base station is 
optimally located at one of the sensor locations. Empirical evidence suggests that 



OO . the result holds in some generality. 

o ; 1 Introduction 

OO 

o 

The recent confluence of low-cost sensing, communication, and computation technologies 
^ i has led to much interest in the development of wireless sensor networks for estimation and 
^ I control. These wireless sensor networks provide for an economical means of extracting 
greater performance and efficiency in a variety of applications, such as manufacturing and 
chemical plant processes [?], indoor climate control (HVAC) [?], environment monitoring 
[?], electrical power distribution [?], and automatic traffic flow [?]. 

There are many advantages that wireless sensor networks hold over their wired coun- 
terparts. For example, the fact that communication is performed wirelessly eliminates 
any physical connection between the various nodes in the network. This enables wireless 
sensor networks to be implemented without any major infrastructure overhauls. This 
also allows for the design of networks that can gracefully accept dynamic changes to its 
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structure, such as the addition or subtraction of nodes, or in the case of mobile sensor 
networks, changing communication topologies. 

While the physical disconnection within the network allows for many practical and per- 
formance related advantages, it also introduces many fundamental issues (e.g. loss/delay 
of information, limits on communication bandwidth, power constraints [?]) that are of 
lesser concern in wired configurations. Issues such as packet loss/delay are crucial in 
estimation and control problems where information flow is assumed to be uninterrupted, 
and much work has been done to examine performance in the presence of these compli- 
cations [?, ?, ?]. 

In this paper, we consider the problem where (noisy) observations of a dynamic process 
are made by multiple sensors at fixed locations. These sensors communicate over a 
lossy wireless network to a base station that computes the optimal state estimate by 
processing the data via a Kalman filter. Since the physical distance between the sensors 
and the base station affects the packet loss probability (which in turn affects the quality 
of the estimates), our goal is to determine the optimal position of the base station. A 
similar problem of locating control logic has been studied in [?]. 

There has been much recent work on decentralized estimation problems where no single 
node is burdened with the majority of the computational tasks [?, ?]. However, imple- 
mentation of these ideas require "smart" sensors that have local computational abilities. 
While these approaches have been shown to possess various performance advantages, 
they are not always realizable. In order to examine more fundamental issues regarding 
the implementation of wireless sensor networks, we restrict our attention to the case 
where stationary sensors merely sense and communicate, and possess no computational 
power. We leave the more general problem for future work. 

The remainder of the paper is organized as follows. In Section [21 we give a brief overview 
of Kalman filtering. Section [3] contains our problem formulation and main results. Some 
discussion regarding more general cases can be found in Section |H 

2 Preliminaries 
2.1 Kalman Filtering 

We begin by giving a brief overview of the canonical state estimation problem. Consider a 
process whose dynamics are modeled by the (possibly time- varying) state space equation 

Xk+i = Axk + Wk, (1) 

where Xk G R" is the state vector at time k, and w is a sequence of independent Gaussian 
random vectors with zero mean and covariance matrix Q. The initial condition xq is 



unknown and modeled as a Gaussian random variable with zero mean and covariance 
matrix Pq. 

We monitor this process via measurements of the form 

Uk = Cxk + Vk, 

where C is possibly time-varying, and w is a sequence of independent Gaussian random 
vectors with zero mean and covariance matrix R. The sequences w and v are assumed 
to be uncorrelated. The minimum variance estimate of the sequence x is then given 
recursively by the Kalman filter [?], 

Kk = APkC*iCPkC* + R)-' 
Xk+i = Axk + Kkivk - Cxk) 
Pk+i = APkA* + Q-KkCPkA 

Here, M* denotes the complex conjugate transpose of a matrix M. 

Due to the increasing availability of high quality, low-cost sensors, we are often con- 
fronted with the situation where we have access to measurements from multiple sensors. 
That is, suppose that we monitor the process ([1]) via sensors, whose measurements 
are represented as 

yl = CiXk + vl 

= CNXk + Vk- 

Here, i;^, . . . , are sequences of independent Gaussian random vectors with zero mean 
and covariance matrices . . . , Rn, respectively. The sequences v^, . . . and w are 
assumed to be uncorrelated. The resulting state estimates are again given by ([2]), where 
we now define 



(2a) 
(2b) 
(2c) 



C = blkcol(Q), j = l,...,N 
R = blkdiag(i?j), j = l,...,N. 

Here, blkcol denotes the vertical stacking of its arguments, and blkdiag denotes a block 
diagonal matrix of its arguments. 



2.2 Kalman Filtering with Packet Losses 

We now discuss the minimum variance state estimation problem in the presence of packet 
losses. Let Q be the power set of {1, . . . , N} and let each subset G O be given the 
standard simple order. Since the Kalman filtering equations ([2]) yield the minimum 



variance state estimate for linear time-varying systems, we obtain the following optimal 
filter for state estimation with packet losses 



and 



where 



Xk+i = Axk + Kkivk - F^^Xk) } for ujk e 
Pfc+i = APkA* + Q- KkF^^PuA* 



Xk+i = Axk 

Pk+i = APkA* + Q 



for uJk 



F^, = blkcol(Cj), j e Uk 
G^^ = blkdiag(i?j), j e u^. 



For a similar treatment, see [?, ?] 



3 Base Station Placement 
3.1 Problem Formulation 

We will model the packet arrival process for each sensor as a sequence of independent 
binomial random variables. More specifically, for i = 1, . . . , A^, let Aj be the probability 
that the observation from sensor i is received. Since the sensor locations are fixed, a 
given base station location gives rise to distances di, . . . , (i^v We will assume that the 
probability of packet arrival decreases as the distance di between a sensor and the base 
station increases. 

Note that the sequences {F^)k and {Gi^)k both depend on whether packets have been 
received or lost, and hence the sequence {Pk)k is itself random. We will then aim to 
minimize (with respect to the Loewner ordering [?], i.e., X y Y if X — Y is positive semi- 
definite), the expected estimation error covariance. That is, we consider the problem 

minE[Pfc+i|Pfc] , 

where the minimization is performed over the set of possible base station locations. If 
we define for each nonempty ^ u & Q 

= n^:'- n (1-^^) 

jGoj je{l,...,N}\uj 

Kk = APkF:{F^PkF: + G^Y\ 



and = for a; = </>, we can write 

As stated, the above optimization problem is a difficult nonlinear programming problem 
that can be attacked by the appropriate solvers. Since it is not feasible for the location 
of the base station to change over time, we will insist that the solution to the above 
problem minimizes E for any Pj.. Consequently, we shall hereafter restrict our 

attention to an illuminating instance of the above problem. 

3.2 The 2 Sensor Case 

In this section, we consider the case where measurements are obtained from 2 homoge- 
neous sensors, Si and S2- Without loss of generality, these sensors are located at = 
and ^2 = 1 on the real line. The goal is to find the best location d e [0, 1] for the base 
station. Note that optimal solution must lie in the interval [0, 1]. The situation ci = 
(or d = 1) corresponds to the base station being physically wired to one of the sensors, 
and wireless communication is performed only by the other sensor. 

Although there may be many factors that influence the packet loss probability, we shall 
consider only the effects due to the physical distance between the communicating nodes. 
Since the sensors have identical broadcasting capabilities, it is then natural to model 
the probability of packet loss in the following manner. Let / : [0, 1] — [0, 1] be convex 
and decreasing, with /(O) = 1. When the base station is located at some d G [0, 1], the 
probability that a packet is received from sensor si is f{d), and the probability that a 
packet is received from sensor S2 is /(I — d). The constraint that /(O) = 1 captures the 
implicit assumption that no packet is lost when a sensor is affixed to the base station 
and communication is not performed wirelessly. 

3.2.1 Vector-valued Measurements, Same Covariances 

We will ffist restrict our attention to the case where 

= C2 = C e R'"^" 
= R^ = Re E™^"*. 

If we model the probability of packet loss as described above, we can examine the 
estimation error covariance as a function of d. For ease of notation, let us define 

Ml = APkC*{CPkC* + R)-^CPkA* 
Ma = APkC*{CPkC* + R/2y^CPkA* 
S = APkA* + Q. 



The estimation error covariance Pk+i then assumes values as follows: 



Pk+l 


Probability 


S-M2 


fid)fil-d) 


S-Mi 


f{d)-2f{d)f{l-d)) + f{l-d) 


S 





Note that in the case where packets from both sensors are lost, the resulting covariance 
corresponds to that of propagating the previous state estimate. The expected error 
covariance can then be computed as 

E[Pfc+i|Pfc] = S- f{d)f{l-d)M^ 

-if id) - 2/(rf)/(l -d) + /(I - d))M2. 



We have the following result. 



Proposition 1 Let Ci = C2 = C e R™^" and = = R e R™x™. Let f : [0, 1] ^ 
[0, 1] be twice differentiable, convex, and decreasing, with /(O) = 1. Then, E [Pfc+i|Pfc] is 
matrix concave in d on [0, 1]. 



Proof Define 



Jk+i{d) = E[Pk+i\Pk] 

= APkA* + Q-f{d)fil-d)M2 

-{f{d)-2f{d)f{l-d) + f{l-d))M,. 

Taking the second derivative yields 

j'Uiid) = r\d)u{i-d)-i)M, 

+/"(!- d)(/(d)-l)Mi 

-2/'(rf)/'(l-rf)(2Mi-M2) 

+ (/" -d) + f{d)f"{l - d)){M, - M2). 

Since / is convex and decreasing with /(O) = 1, we have f'{x) < and f"{x) > for 
X G [0,1]. Consequently, J"{d) ^ 0, and hence E[Pfc+i|Pfc] is matrix concave on [0,1] 
(see [?], p. 110). ■ 



Corollary 1 Let Ci = C2 = C E E™^" and R^ = R^ = R e R™^'". Let f : [0, 1] ^ 

[0, 1] be twice differentiable, convex, and decreasing, with /(O) = 1. Then, the base 
station is optimally located at one of the sensor positions. 



3.2.2 Scalar-valued Measurements, Same Covariances 



A similar result may be obtained in the case where Ci = C2 = C G R^^" and the noise 
covariances -Ri and R2 are not necessarily equal. We will need the following preliminary 
result. 

Lemma 1 Let Ci = C2 = C e E^^" and define 
T = CPkC* 

n * r n —1 

Mo = AP Ft^' A rjPA* 
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'T + Ri T 
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'C 


T T + R2 
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Ml = 
M2 = 

Then, Mo - Mi - M2 h 0. 
Proof Let us define 



APC*{T + Ri)-^CPA* 
APC*{T + R2y^CPA*. 
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can be written as 
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and U 
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where the second equality follows from the matrix inversion lemma [?]. Now, let 

(1 - 7-^/5)' 



H 



a 



(1- 7-^/3)^ -a-Ha-/3V^) 



The desired result follows from noticing that the denominator is positive and the nu- 
merator satisfies 



= ( 2 

< 0. 



CPuC* 



CPkC* 



CPkC* + i?2 CPkC* + Ri 



Proposition 2 Let Ci = C2 = C G R^^"" and suppose that Ri and R2 are not neces- 
sarily equal. Let f : [0, 1] [0, 1] be twice differentiable, convex, and decreasing, with 
/(O) = 1. Then, E [Pfc+i|Pfc] is matrix concave in d on [0, 1]. 



Proof Define 

T 

Mo 

Ml 
M2 

We can then write 

Jk+i{d) 
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c 



PA* 



APC*{T + RiY^CPA* 
APC*{T + R2)-^CPA*. 



E[Pk+i\Pk] 
APkA* + Q 

-f\d)f{l - d)M, - f{d){l - /(I - rf))Mi 

-j\l-d){l-f\d))M2 



Taking the second derivative yields 

JUM = nd){f{i-d)-i)M, 

+/"(!- rf)(/(d)-l)M2 
-2/'(rf)/'(l-rf)(Mi + M2-Mo) 
+r(rf)/(l-rf)(M2-Mo) 
+/(rf)r(l-rf)(Mi-Mo). 

Since / is convex and decreasing with /(O) = 1, we have f'{x) < and f"{x) > for 
X e [0, 1]. From Lemma [1] we have Mi + M2 - Mq ^ 0. As a resuh, J"{d) ^ 0, so 
E [Pfc+il-Pfc] is matrix concave on [0, 1] (see [?], p. 110). ■ 



Corollary 2 Let Ci = C2 = C E R^^*^ and suppose that Ri and R2 are not necessarily 
equal. Let f : [0, 1] [0, 1] be twice differentiable, convex, and decreasing, with /(O) = 1. 
Then, the base station is optimally located at one of the sensor positions. If Ri >z R2, 
then the base station is optimally located at the second sensor. 

3.2.3 An Example 

We now offer an illustrative example to complement our main results. 

Suppose that there are two sensors, located respectively at the points and 1 on the 
real line. Let us model our packet arrival probability function as 

f{d) = e-". 

Note that / is convex and decreasing on [0, 1] with /(O) = 1. The process being mon- 
itored is described by the state space equations in ([T]), where all the eigenvalues of 
A G R^^^ lie in the unit circle and C G R^^^. The trace of the estimation error covari- 
ance matrices is plotted in Figure [U for the case where the base station is situated at 
the point and the case where it is situated at the point 0.5. 
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Figure 1: Covariances for = (dashed) and d = 0.5 (solid). 

4 General Cases 

We now offer some commentary on extending our main results. 

Although the previous results were confined to the case Ci = C2, simulations suggest 
that they also hold in more generality. Consider the case where N sensors located 



colinearly are available to make observations. In general, we may have different values 
for Ci, . . . , Cat and different values for Ri, . . . , R^. In this case, empirical evidence 
supports the extension of our main result in the sense that the expected error covariance 
is piecewise concave over the intervals defined by the location of the sensors. Figure [2] 
illustrates this situation for the case N = 3. 
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Figure 2: Piecewise concavity for multiple colinear sensors. 




In the case where the convex hull of the sensor locations has a nonempty interior in R^, 
it can be shown that the optimal base station location may lie at a point other than the 
sensor locations. Consider 3 sensors, located at the vertices of an equilateral triangle 
with side length 0.5. For = /, f{d) = e~^, Ci = C2 = C3 = 0.6, and unit noise 
covariances, the expected error covariance (as a function of base station location) has a 
local minimum at the centroid, as illustrated in Figure [31 
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Figure 3: Local minimum in the interior of the convex hull. 



5 Conclusion 



We considered the problem of positioning a central computational node within a wireless 
sensor network for the purpose of state estimation. In the presence of potential packet 
losses, the optimal estimate is given by a time-varying Kalman filter. We then strived 
to determine the optimal location of the base station so as to minimize the expected 
conditional estimation error covariance. It was shown for the case of a 2-sensor config- 
uration that under a certain class of packet loss probability models, the base station is 
optimally located at one of the sensor positions. 



